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Introduction. The integration of reliability and optimization 
concepts seeks to design structures that should be both 
economic and reliable. This model is called Reliability-Based 
Design Optimization (RBDO). In fact, the coupling between 
the mechanical modelling, the reliability analyses and the 
optimization methods leads to very high computational cost 
and weak convergence stability. 

Materials and Methods. Several methods have been developed 
to overcome these difficulties. The methods called Reliability 
Index Approach (RIA) and Performance Measure Approach 
(PMA) are two alternative methods. RIA describes the 
probabilistic constraint as a reliability index while PMA was 
proposed by converting the probability measure to a 
performance measure. An Optimum Safety Factor (OSF) 
method is proposed to compute safety factors satisfying a 
required reliability level without demanding additional 
computing cost for the reliability evaluation. The OSF 
equations are formulated considering RIA and PMA and 
extended to multiple failure case. 

Research Results. Several linear and nonlinear distribution 
laws are applied to composite yarns studies and then extended 
to multiple failure modes. It has been shown that the idea of 
the OSF method is to avoid the reliability constraint evaluation 
with a particular optimization process. 

Discussion and Conclusions. The simplified implementation 
framework of the OSF strategy consists of decoupling the 
optimization and the reliability analyses. It provides designers 
with efficient solutions that should be economic satisfying a 


required reliability level. It is demonstrated that the RBDO 
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Beedenue. Vunrerpauwa KkonlenuMii HagexKHOCTH U 
ONTHMUH3alIM HalipaBsieHa Ha pa3spaOoTKy CTpyKTyp, KOTOpbIe 
JOJDKHBI ObITb JKOHOMHYHEIMM WM HajexKHbIMH. DTa MOJesIb 
Ha3biBaeTca OMNTHMH3al[MA TpOeKTHPpOBaHHA Ha OCHOBEe 
(RBDO). Me@KILy 


Me€XaHH4eCCKHM MOJICJIMPOBAaHHeM, aHaJIH3OM Hajl@oxKHOCTH HU 


HayleoKHOCTU @Maxtuyecku CBA3b 
MeTOJaMH ONTHMH3aUuHH NpHBOAHT K O4CHbB BbICOKHM 


BEIYHCIMTCIbHBIM 3aTpaTamM wu ciladol cTadMIbHOCTH 
KOHBepreHuHu. 

Mamepuaavi u memoooi. Ana ipeoqonenuad 9THX TpyqHOcTeli 
Obuio paspaboTaHO HecKoNbKO MeToyoB. MetTogb nog 
Ha3sBaHHem «MuHyekc mMoKa3aTeia HafexKHocTH» (RIA) u 
«Metoq oeHKH o92pdextupHoctu» (PMA) - 9To gpa 
ayIbTepHATHBHBIX MeToya. RIA onmMcbrBaeT BeposxTHOCTHOe 
orpaHw4eHve Kak HHJeKC HaexKHOCTH, B TO Bpema Kak PMA 
OBL TpeqoxwKeH Kak HyTb mpeoOpa3s0BaHHA BepOATHOCTHOI 
Mepbl 23dpdextuBuoctu. IIpeqnoxeH MeTOX ONTHMasIbHOrTO 
Kooppuynenta Oe3omacHoctm (OSF) gma  pacueta 
Koo) PHUMeEHTOB 
TpeOyeMoMy YpoOBHIO Hajle@xKHOCTH, 6e3 OMOJHHTebHBIX 


BbIGHCJINTCJIBHbBIX 3aTpaT JUJLA OLCHKU 


Oe30macHocTH, YHOBJICTBOPALOLIAX 
HaJl@KHOCTH. 
Ypasuenua OSF ccibopmynupopansr c yyetomM RIA u PMA u 
PacluMpeHbl JO cily4ad C HECKOJIBKAMH OTKa3aMH. 

Pe3ynomamot ucciedoganua. HeckombKo JIMHeMHbIX 4H 
HeJIMHCHHBIX 3AKOHOB paciipeyeyieHHa pMMeHAIOTCA K 
MCC€HOBaHHAM KOMIIO3HIMOHHBIX HUTCH, a  3aTeM 
PactipocTpaHsKOTCA Ha MHO2KECTBCHHbIC Pe2KHMbI OTKA3a. 
Baio mo0Ka3aHoO, YTO Hed MeToya OSF 3aksrouaeTca B TOM, 
YTOOBI H30e%KATb OLCHKH OrpaHW4eHHA HalexKHOCTH Cc 
TIOMOLMIbIO KOHKpeTHOLO Mpoljecca ONTHMUH3alHH. 

O6cyorcoenue u eoigodbl. YIpouleHHas CTpyKTypa peasu3alyun 
ctpateruH OSF 3axmouaeTca B pa3qelleHHH aHamnH3a 
OnTHMv3alMH «6M CHayjexKHOCTH. Ona  lpeocraBiaeT 
TIpOCKTHpOBIHKaM  9(PpdPeKTHBHEIe KOTOpBIe 
OJDKHEI JKOHOMHYECKH YOBICTBOPATh TpeOyeMbIit ypOBeHb 


pelieHua, 
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compared to OSF has several advantages: small number of 
optimization variables, good convergence stability, small 


computing time, satisfaction of the required reliability levels. 
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1. Introduction 

When Deterministic Design Optimization (DDO) methods are used, deterministic optimum designs are usually 
pushed to the design constraint boundary, leaving little or no room for tolerances (or uncertainties) in design, 
manufacture, and operating processes. So deterministic optimum designs obtained without consideration of 
uncertainties could lead to unreliable designs, therefore calling for Reliability-Based Design Optimization (RBDO). It is 
the objective of RBDO to design structures that should be both economic and reliable. However, the coupling between 
the mechanical modeling, the reliability analyses and the optimization methods leads to very high computational cost 
and weak convergence stability. To overcome these difficulties, two points of view have been considered. From a 
reliability view point, RBDO involves the evaluation of probabilistic constraints, which can be executed in two different 
ways: either using the Reliability Index Approach (RIA) or the Performance Measure Approach (PMA) [1]. The major 
difficulty lies in the evaluation of the probabilistic constraints, which is prohibitively expensive and even diverges for 
many applications. However, from an optimization view point, a hybrid method based on simultaneous solution of the 
reliability and the optimization problem has successfully reduced the computational time problem. However, the hybrid 
RBDO problem is more complex than that of deterministic design and may not lead to local optima. To overcome both 
drawbacks, an Optimum Safety Factor (OSF) method has been proposed to compute safety factors satisfying a required 
reliability level without demanding additional computing cost for the reliability evaluation. The efficiency of the OSF 
method has been demonstrated relative to the hybrid one for a linear distribution law [2]. In this work, the OSF is 
reformulated considering RIA and PMA. Next, an extension to multiple failure modes is carried out. Finally, an 
application on composite yarns is carried out for linear and nonlinear distribution laws and multiple failure modes 
(scenarios). 

2. Reliability-Based Design Optimization 
The reliability-based design optimization problem is performed by nesting two sub-problems [3, 4]: 


1: Optimization problem: The objective is to minimize an objective function f (2) subject to deterministic 
constraints g, (x) <0 and a required reliability level B(x,u) > B, as follows: 

min : f (x) subject to g, (x) <0 and B(x,u) =f, (1) 

2: Reliability problem: The objective is to find the minimum distance between the origin of the normalized 


space and the MPP (Most Probable failure Point or design point) on the limit state function (Fig. 1). The problem can be 
written as follows: 





n 
ye subject to H(x,u) <0 (2) 
1 
For more details about reliability analysis, readers can see [5]. These two sub-problems are carried in two 
different spaces: physical and normalized spaces (Fig. 1). 
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Fig. 1. Physical and normalized spaces [3, 4] 


These nested problems require a high computing time, however, the Optimum Safety Factor (OSF) is used to solve 
these numerical difficulties. 

3. OSF developments 

In the RIA and PMA, the major difficulty lies in the evaluation of the probabilistic (reliability) constraints; this 
is prohibitively expensive and even diverges for many applications [6]. RIA describes the probabilistic constraint as a 
reliability index while PMA was proposed by converting the probability measure to a performance measure. It has been 
found in the literature that PMA is more efficient and stable than RIA in the RBDO process [1]. Furthermore, PMA 
converges to different optimum solutions when starting from different initial designs [7]. In this section, it is shown that 
the OSF formulations can be deducted from both approaches (RIA & PMA). 

3.1 OSF considering RIA and PMA 

3.1.1 OSF based on RIA 

The probabilistic constraint in RIA can be evaluated by solving the first-order reliability analysis to calculate 
the design point P*, which is formulated as an optimization problem with an equality constraint H (u) = 0, in the 
normalized space: 








min : d(u) 
7 (3a) 
subject to: H(u)=0 
and can be also written as: 
: 2 
min > d“(u 
u si (3b) 
subject to: H(u)=0 
The Lagrangian function for the problem (3b) can be written as 
L(u,A)=d?(u)+2-H(u) (4) 
The optimality conditions for the Lagrangian function are: 
2 
ae oh a0 , d=h..n (Sa) 
Ou; Ou; Ou; 
OL 
—-=H(u)=-0 5b 
au) (5b) 
Using the expression for the square distance d’ in the equation (3b), we get: 
ype, i=l,...,.n (6) 
2 Ou; 


3.1.2 OSF based on PMA 
The evaluation of the probabilistic constraint in PMA requires an inverse reliability analysis, which corresponds 


to the inverse problem of the reliability analysis to calculate the design point P*, which is formulated as an optimization 
problem with an equality constraint d : (u) = B? , in normalized space: 
min : Hu) 
we (7a) 
subject to: d(u) =f, 


and can be also written as: 
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min : H(u) 
i (7b) 
subject to: d” (u) = B? 
The Lagrangian function for the problem (7b) can be written as 
L(u,2,s)=H(u)+A-|d?(u)-B? | (8) 
The optimality conditions for the Lagrangian function are: 
2 
ay 5 ies (9a) 
Ou; Ou; Ou; 
OL 
—=d*(u)=0 9b 
On ( ) a 
Using the expression for the square distance d’ in the equation (7b), we get: 
uj et ; i=l,....n (10) 
2r Ou; 


According to equations 6 and 10, the resulting normalized vector u can be written in term of the same 
derivative 0H /0u; . 


3.2 OSF for component RBDO 
Let us consider now the case of n=2 normalized variables i = /, 2 (see Fig. 2). The tangent ofais given 


by: tana =u, /u,. Using (6) and (10), we get the same formulations as follows: 
oH 


Ouy 
an H ( ) 


du 
The formulation (11) shows that there is no difference between RIA and PMA for the following OSF developments. 
The problems (3) and (7) give us the reliability index 8 as the minimum distance between the limit state surface and the 


origin [8]. This means that the resulting reliability index may be lower or higher than the target reliability index B,. As 


we wish to satisfy a required target reliability level for the optimization problem, we can write 


Bp= ou (12) 
i=l 


Thus, the general expression for the normalized variable u; when using RIA or PMA, is 


, i=h.n (13) 








The calculation of the normalized gradient OH’ /0u is not directly accessible because the mechanical analysis is 
carried out in the physical space rather than in the standard space [9, 10]. However, using theory of statistics, we can 
derive the following expression from which the computation of the normalized gradient can be carried out by applying 
the chain rule on the physical gradient 0G / Oy : 


dH _ G O%'(x,u) 
Ou; Oy, = Ou; 





,i=l,...7, k=1,...,.K (14) 


where 7'(x,u) is the inverted probabilistic transformation function. We find that the normalized gradient can be 


expressed as 
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0H {0G 
Ou; Oy; 





,i=L.n (15) 


To summarize, we satisfy the required target reliability index as follows: 


f=L.n (16) 





Thus, we compute the optimum safety factors subject to the optimum values of the normalized variables u, . According 


to (13) and (15), it has been demonstrated that the optimum values of u; can be written in the following form: 


i=l,...n (17) 








where the sign of + depends on the sign of the derivative, i.e.: 


OF Gu 30 and Sw 20- 4 than (18) 


i i 





Fig. 2. Design point modeling in a double normalized variable space 


The efficiency of the OSF is to transfer the reliability problem into an analytical formulation (17). This 
formulation has been extended to several distribution laws (normal, lognormal, uniform Weibull, and Gumbel laws) for 
single failure mode cases considering only RIA [11]. However, in this work, both RIA and PMA are used to show the 
efficiency of the OSF formulation, and an extension to multiple failure case is carried out in the next section. 

3.3 OSF for system RBDO 

The system reliability problem can be written as: 








=min d(u;) = Jup +ug +...+u2 st: H;(u,u2,...um_) <0 (19) 


Bete i 


where d(u;) is the minimum distance between the Most Probable failure Point (MPP) and the optimal solution 
in the normalized space (Fig. 1). And H; (u; ) <0 represent the different failure modes. Using similar way, the 


analytical formulation using OSF for a system RBDO can be written as follows: 
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(20) 





The formulation (20) provides different optimum values of the normalized variables at the failure point and 
taking into account several failure modes. 

4. Numerical Application on composite yarns 

Composite materials have been used in structures for centuries. In recent times, composite parts have been 
used extensively in aircraft structures, automobiles, sporting goods, and many consumer products. Composite materials 
are those containing more than one bonded material, each with different structural properties. The main advantage of 
composite materials is the potential for a high ratio of stiffness to weight. To solve the composite yarns structures (or 
systems), a stress unilateral study is required. 





M2 








M1 











b) Section of the studied rope 


a) Ropes of composite yarns 


Fig. 3. Material structure 


4.1 Problem description 
For tri-material structure illustrated in Fig. 3, the number of yarns n, , the section S;, the Young's modulus E£,, 


the Poisson's ratio v,and material density p; are given according to each material as shown in Table 1. 














Table 1 
Material properties of different used yarns 
Material number n; S, (m’) E, (Pa) V; p; (Kg/m*) 
M1 5 10e-8 5.86 0.4 4.0 
M2 15 10e-8 8e6 0.3 5.0 
M3 10 10e-8 9e6 0.2 6.0 


























4.1 Single yarn behavior 

A stress unilateral structure is not capable of transmitting compression, i.e. negative stress: all its internal efforts 
are traction. The formulation of this problem is complex and may lead to non-existence of solutions or several 
approximate solutions. Some methods have been developed to determine a solution among all resulting configurations 
[12]. 





= o lll 
: p=p,(l+e) 


a) Initial case of a yarn and b) Stress case when applying external forces 


Fig. 4. The tension on the yarn 
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The behaviour of an elastic file is not defined by the single Hooke’s law; even in the frame of the linear elasticity, we 
have 


jot and T>0 (21) 
K 


where K >0 is the constant of elasticity and T is the tension which defines the internal efforts of the studied yarn. The 
total section is given by 
n3 


n n 
a > sy +s, +> S34 (22) 
i=l j=l k=l 


and the total tension is given by 
nQ n3 


ii Sh, ++ > Be (23) 
a] ‘=| k=l 


with 
Sti So j S34 
T.=—T,T,,=—T and 7,, =—=T 
Li S 2,J S 3,k S 

Fig. 4a shows the initial case of a yarn in which the volume mass =p). When applying an external force and 
considering the gravity force (see Fig. 4b), the equilibrium vector equation can be written as 

dT 

—+pg¢=0, (24) 

da 
with p=p,(.+6). 
For our studied rope, there is a non gravity force ppg = 0 . So, we get 

T(L) = fox (25) 

This way, the tension is a proportional value of the external force as follows 





Rc Shi de S97 dlc S34 be 
Li = Sex 2,7 — S See an 3,k — S Tee ( ) 
Here, the failure conditions are written as follows: 
a1 s £1 max > £2 =< £2 max and £3. =< ©3 max (27) 
In order to compute the maximum strain of each material, we have 
Tj Th; Ts, 
eo dso ond oe Se (28) 
li K 25 3,k K 
Li 2, 3k 


4.3 Procedures 


DDO procedure: 
The optimization problem is to minimize the mass subject to the allowable strain ¢,, and taking a global 


safety factor S, . This optimization problem must be followed by the reliability analysis as follows: 
1- Optimization problem 
min mass(x) subject to (x) <&,, = €1im/ Sy (28) 


2- Reliability analysis 


subject to e(u) < €;,, (29) 





RBDO procedure 
Using OSF method, the RBDO procedure contains three main steps: 


1- The first step is to obtain the design point (MPP). It is the objective to minimize the volume subject to the design 
constraints without consideration of the safety factors. This way, the optimization problem for a single failure mode is 
simply written as: 


min mass (y) 
subject to e(y) < &jin Co 


and for a multiple (double) failure mode as follows: 


min mass (y) 
subject to e(y) < &jjm (31) 
T(y) s Tim 
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2- The second step is to compute the optimum safety factors using linear and nonlinear distributions [9] when the 
number of the deterministic variables is equal to that of the random ones. During the optimization process, we obtain 
the sensitivity values of the limit state with respect to all variables. 
3- The third step is to calculate the optimum solution. This encompasses inclusion of the resulting values of the safety 
factors into the design variables in order to evaluate the optimum solution. 

4.4 Numerical results 

4.4.1 Component RBDO 

According to Table 2, the DDO procedure leads to a high reliability index, while the RBDO one satisfies the 
required reliability level for the linear and nonlinear distribution laws. The resulting reliability index is: B =3 that 


corresponds to a probability of failure (Py ~ 0.1% ). 






































Table 2 
DDO and RBDO results for a single failure mode 
Desi Optimum Solution 
Parameters as DDO RBDO 
Normal Normal Lognormal Uniform Weibull 
Sl 0.2900 0.4400 0.3415 0.3387 0.3413 0.3455 
S2 0.4750 0.8500 0.5774 0.5698 0.5655 0.5922 
S3 0.5550 1.2300 0.6846 0.6737 0.6632 0.7072 
Cian 0.0095 0.0048 0.0078 0.0079 0.0080 0.0076 
Mass 0.034 0.0669 0.042 0.041 0.041 0.043 
B -- 7.8231 3.00 3.00 3.00 3.00 
Sy -- -- 0.849274 0.856135 0.849613 0.839352 
Spo -- -- 0.822641 0.833690 0.839981 0.802099 
S93 -- -- 0.810724 0.823839 0.836908 0.78483 1 





























4.4.2 System RBDO 

It is easy to study structures under the most critical mode but it may provide inaccurate results for structural 
reliability analysis. When dealing with several failure modes, the sensitivity study of each failure mode with respect to 
all parameters may lead to strange results for the role of certain parameters. The required reliability for both failure 
modes (tension and strain) is considered to be: B, =3 that leads to a system reliability index B.s,.,.,, 2B, . The system 


reliability index should be bigger or equal to the minimum required reliability index of all failure 
modes B5,.em 2 min(B;). Table 3 shows Linear and nonlinear RBDO results for the same target reliability index for 


both failure modes. 






































Table 3 
Linear and nonlinear RBDO results for multiple failure modes 
Panes DESH Optimum pobition 
Point Normal Lognormal Uniform Weibull 
SI 0.2900 0.3432 0.3402 0.3421 0.3479 
S2 0.4750 0.5774 0.5698 0.5655 0.5922 
S3 0.5550 0.6850 0.6740 0.6632 0.7078 
€ 0.0095 0.0078 0.0079 0.0080 0.0076 
T 66335 58956 59526 59909 57830 
Mass 0.0343 0.0418 0.0413 0.0409 0.0430 
B ation 0.0000 3.0242 3.0242 3.0241 3.0183 
Sy 0.00 0.845081 0.852562 0.847811 0.833615 
Sy 0.00 0.822641 0.833690 0.839981 0.802099 
S93 0.00 0.810263 0.823460 0.836803 0.784156 
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6. Conclusions 

In this paper, the OSF method is shown as a distinctive tool for RBDO problems. This method is based on the 
sensitivity of the limit-state function with object of determining the role of each studied parameter relative to the failure 
mode (or modes). First of all, it has been shown that the idea of the OSF method is to avoid the reliability constraint 
evaluation with a particular optimization process. In addition to its simplified implementation framework to completely 
decouple the optimization and the reliability analyses, it provides designers with efficient solutions that should be 
economic satisfying a required reliability level. The developed equation of OSF basing on both of RIA and PMA leads 
to the same formulations. The OSF procedure needs only single optimization process for the design point without 
additional computing time because it has a single variable vector. Finally, the developments based on the reliability 
view point are less efficient than those based on the optimization view point because the second provides us with 
reliability-based optimum designs without additional computing cost for probabilistic (reliability) constraints and can 
lead to global optima. It is shown that the RBDO compared to OSF has several advantages: small number of 
optimization variables, good convergence stability, small computing time, satisfaction of the required reliability levels. 
For composite structures, the RBDO problems are more difficult than for simple structures because we deal with several 


related limit state functions. 
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